Cosmography and cosmic acceleration 
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We investigate the prospects for determining the accelerating history of the Universe from upcom- 
ing measurements of the expansion rate H(z). In our analyses, we use Monte Carlo simulations based 
on wCDM models to generate samples with different characteristics and calculate the evolution of 
the deceleration parameter q(z). We show that a cosmographic (and, therefore, model-independent) 
evidence for cosmic acceleration (q(z < zt) < 0, where Zt is the transition redshift) will only be 
possible with an accuracy in H(z) data greater than the expected in current planned surveys. A 
brief discussion about the prospects for reconstructing the dark energy equation of state from the 
parameters H(z) and q(z) is also included. 

PACS numbers: 98.80.-k, 95.36.+X, 98.80.Es 



I. INTRODUCTION 



The determination of cosmographic parameters, such 
as Hq and qo, has a long and interesting history in Cos- 
mology (see, e.g., In particular, the evolution of 
such parameters provides a unique and direct method to 
map the expansion history of the Universe in a model- 
independent way. Since all evidence we have so far for 
the current cosmic acceleration are indirect 0, extract- 
ing the evolution of these two parameters from future 
redshift surveys constitutes one of the major challenges 
in observational cosmology 1 . 

In this paper, we investigate how well cosmography 
may provide a model-independent way to check the real- 
ity of cosmic acceleration. Specifically, we study the evo- 
lution of the deceleration parameter from the H{z) data 
which are to become available by some planned projects. 
To this end, we performed a Monte Carlo (MC) simula- 
tion based on wCDM models (of which the ACDM model 
is a special case with w = p/p = —1, where p stands for 
the dark energy pressure and p for its energy density). 
We generated three samples of H(z) with increasing ac- 
curacy and derived q(z) from numerical differentiation. 
We show that a comographic evidence for cosmic accel- 
eration will be possible only with a great (maybe out of 
the perspective of current planned surveys) accuracy in 
H(z) data. It is worth mentioning that the determina- 
tion of q{z) will also allow us to estimate the transition 
redshift, z tl which corresponds to the epoch when the 
Hubble expansion switched from a decelerating to an ac- 
celerating phase. Since q(z) oc w(z)H(z) 2 / f(z,w), the 
same technique can also be used to map the recent past 



evolution of w(z) in a given model of dark energy. 



II. DIRECT OBSERVATION OF H(z) 

Recently, it has been shown that Luminous Red Galax- 
ies (LRG's) can provide us with direct measurements of 
the expansion rate H(z) [3] (see also Q for a recent re- 
view on H{z) measurements from different techniques) 2 . 
This can be done by calculating the derivative of cosmic 
time with respect to redshift 
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from measurements of age difference between two pas- 
sively evolving galaxies at different z. Simon, Verde and 
Jimenez [7j (SVJ) have demonstrated the feasibility of di- 
rectly measuring H(z) using this differential age method 
and have provided 9 determinations of the expansion rate 
at z ^ 0. More recently, Stern et al. [8| have updated 
SVJ sample to 11 estimates of H(z) lying in the red- 
shift interval 0.1 < z < 1.75. New age-redshift datasets 
for different galaxy velocity dispersion groups (7 LRG's 
sample) have also been made available by the SDSS col- 
laboration n. 

Ref. Q have also pointed out that in the near future, 
the Atacama Cosmology Telescope (ACT) 3 is expected 
to provide observations of over 500 galaxy cluster up to 
z < 1.5. This, together with spectra to be acquired in 
other telescopes in Chile and Southern African Large 
Telescope (SALT) in South Africa, will provide a sam- 
ple of more than 2000 passively evolving galaxies in the 
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1 Due to the multiple integrals that relate cosmological parameters 
to cosmological distances, direct determinations of H (z) may also 
reduce the so-called smearing effect that makes constraining the 
dark energy equation of state (EoS) w extremely difficult [jj]. 



2 Direct measurements of H(z) at different redshifts will also be 
possible through measurements of the line-of-sight or radial com- 
ponent of baryonic acoustic oscillations (BAO) from large red- 
shift surveys with redshift precision of the order of 0.003(1 + z) 
(see, e.g., Q). 

3 http:/ /www. physics. princeton.edu/act/ 
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FIG. 1: a) A MC realization of 1000 simulated values of the Hubble parameter with 15% accuracy based on Ref. 0. b) 
Simulated values of H(z) shown in the previous panel averaged over bins with width Az = 0.1. Squares represent current H{z) 
measurements of Ref. [g . c) Similar to Panel (a) for the 3% projection made in Ref. [l(J . 



redshift range < z < 1.5. From these observations, it 
will be possible to determine ~ 1000 values of the Hub- 
ble parameter at a 15% accuracy level if the ages of the 
galaxies are estimated with 10% error. 

Following a similar approach, Crawford et al. [Io| ex- 
amined the observational requirements to estimate H (z) 
to a given precision. In their simulation, observations 
of LRG's are made at two redshifts, namely, z = 0.32 
and z = 0.51 (average 0.42) and they supposed that un- 
certainty on ages of individuals galaxies lies in the range 
0.05 — 2 Gyr. They estimated that the uncertainty on the 
mean ages will be 0.10, 0.05 and 0.03 Gyr and that Hub- 
ble parameter can be measured with a precision of 10%, 
5% and 3%. A study was also made on the constraints 
upon the exposure time per galaxy and number of galax- 
ies observed. They concluded that with a total time of 
17, 72 and 184 hours, observing 80, 327 and 840 galax- 
ies on the Southern African Large Telescope, will make 
possible to recover H{z) to 10%, 5% and 3% accuracy, 
respectively fioj (private communication). 



III. SIMULATED DATASETS AND q(z) 

In our analyses, we perform Monte Carlo simulations 
which provide us with samples of H(z) based on the 
flat ACDM model. For each "measurement", we as- 
sume a Gaussian distribution of H(z) centered at the 
value predicted by a flat wCDM with w = — 1 with a 
standard deviation corresponding to the percentage ac- 
curacy predicted in future experiments. We shall exam- 
ine two cases based on the projections made in Refs. Q 
and [l(j. For the former case, we simulated 1000 data 
points for H(z), uniformly distributed between z = 
and z = 1.5 with 15% precision. In the latter case, we 
simulated 15 measurements of the Hubble parameter be- 
tween redshifts and 1.5 and equally spaced (Az = 0.1). 
We take the errors estimated in Ref. [l(| when recov- 
ering H(z) from LRG's ages, namely 10%, 5% and 3%. 
Note that in the last case (the one explored throughout 



this paper), one needs 8400 LRG's in the whole interval 
z = 0.1 — z = 1.5. In all our simulations, we have adopted 
Hq = 74.2 km/s/Mpc, which corresponds to the central 
value given in Ref. [ll| based on differential measure- 
ments of Cepheids variable observations and Vt m = 0.27, 
as given by current CMB measurements [12|. 

Figure la shows one realization of 1000 simulated val- 
ues of H(z) with 15% accuracy according to Ref. Q for 
z in the range (0.0 — 1.5). The binned data points of 
this realization (with a bin width Az = 0.1) are shown 
in Figure lb. The error bars correspond to the standard 
deviation around the mean of individual values of H(z) 
within each bin, which contains approximately 67 galax- 
ies. For comparison, the observed values given in [8] arc 
also shown. In Figure lc we show the results of one MC 
realization taking the errors estimated in Ref. [l(| with 
the resulting H(z) data points at 3% accuracy. 

From the above H(z) simulated data, one may derive 
evolution of the deceleration parameter q(z), defined as 
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In order to do that one needs to compute the derivative 
of the H(z) with respect to the redshift. However, while 
integration tends to smooth out data fluctuations, nu- 
merical differentiation tends to magnify errors and the 
scatter of the points. Thus, depending on the degree of 
scatter of H(z) measurements, the gradient estimation 
may be useless. For equally spaced points, the deriva- 
tive is calculated using finite difference approximation, 
that is, dH(zi)/dt ~ [H(z i+ i) — H(zi-i)]/2Az, where 
Az = Zi — Zi-\. We then substitute this into @ and 
calculate q(zi) and the associated uncertainty by using 
the standard error propagation method. 

For our first simulation, shown in Figures la and lb, 
we use the binned points since the large scatter of the 
raw data makes them useless. The resulting values of 
the deceleration parameter from the MC realization of 
Figure lb are shown in Figure 2a. Although the points 
lie near the predicted curve of the ACDM model (black 
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FIG. 2: The evolution of q(z) derived from the H(z) points displayed in Figure 1. Panel a corresponds to the 1000 binned 
values of Fig. lb whereas Panel b to the 15 H(z) data points at 3% accuracy shown in Fig. lb. c) The difference between the 
simulated and the background model values of q{z) for 1000 MC realization and 15% accuracy in H(z) observations, d) The 
same as in Panel c for the 15 H(z) data points at 3% accuracy shown in Fig. lb. 




FIG. 3: a) The results of our simulations following [lfj] with a precision of 1% in the H(z) observations, b) The evolution of 
q(z) for the MC realization shown in Panel a. With this precision in the H(z) measurements, a cosmographic detection of of 
the current cosmic acceleration can be obtained, as can also be seen from Panel c. 



curve), we cannot expect a definitive evidence for a re- 
cent cosmic acceleration from these simulated H(z) sam- 
ple with 15% accuracy. Note that the same is also true 
even when the evolution of q{z) is derived from the 15 
data points with 3% accuracy following [l(| (Fig- 2b). 
In this case, although the error of individual points is 
smaller, the determination of the deceleration parameter 
is less precise with considerable scattering of points and 
error bars allowing both a decelerating and accelerating 
universe today. 

In order to quantify the above results, we also per- 
formed 1000 MC realizations to calculate the mean of the 
absolute values of the difference between the calculated 
and the background model values of Qacdm (z% ) i that is, 

1 N 

< Aq >= — ^ \lsim{Zi) - <7ACDm(^)|- (3) 

i=l 

We plot this quantity for both cases discussed above in 
Figures 2c (1000 binned points) and 2d (15 points), with 
the error bars representing the standard deviation from 
the mean. As we see, < Aq > has a weak dependence 
on redshift for the first case and in the interval 0.1 < 
z < 1.4 it lies in the range ~ 0.1 — 0.25, although the 



uncertainty is relatively large. From Figure 2d we see 
that the departure of the calculated q S im{zi) from the 
expected value given by the model is larger than in former 
case, with < Aq > varying from 0.2 to 0.4. This clearly 
reflects the large scattering of points shown in Figure 2b. 

In view of the above results, we also speculated about 
how precise should future observations be in order to pro- 
vide a clear evidence for cosmic acceleration from cos- 
mography. We did that for the method proposed in [l(| 
and found that it will give very good results if the accu- 
racy is improved to as much as 1%. This is illustrated 
in Figure 3. In this case, the calculated value of q S im{zi) 
is very close to the ACDM model predictions in almost 
the entire range of z and the deviation from the expected 
value is as low as 0.07 and always < 0.15. Naturally, the 
observational requirements for this case would be much 
more stringent than the previous ones, with the need of 
thousand galaxies and more observational time as well. 



IV. EQUATION-OF-STATE PARAMETER 

By considering a dark-energy dominated universe, a 
step further in the above discussion concerns the recon- 
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FIG. 4: Evolution of w with redshift for 1% accuracy in H(z) 
observations following the method of Ref. [jjj. Although the 
background model is fairly recovered, the uncertainties in- 
crease considerably at the upper end of the redshift interval 
considered. 

struction of its equation-of-state parameter w from mea- 
surements of H(z). In order to be consistent with the 
simulations presented earlier, we consider accelerating 
wCDM models with w lying in the interval [-1/3, -3/2]. 
In this context, Eq. @ can be rewritten as 

[ ' 3n°(l + *)3(i+«0 \Ho) ' ( ) 

where is the current fractional contribution of dark 
energy to the critical density. Note that such procedure 
does not involve several integrations of w(z), as in cos- 
mological tests based on age or distance measurements 
and, therefore, may reduce the smearing effect on w{z) 
determinations discussed earlier (see also [4j for a differ- 
ent approach). 

Figure 4 shows the evolution of w with redshift cal- 
culated from the simulated measurements of H(z) using 
the method proposed of Ref. with 1% accuracy. Al- 
though the input model [w = — 1) is fairly recovered, we 
note that the estimated errors of the dark energy EoS 
(Sw) increase considerably at the upper end of the red- 
shift range. This can be more easily understood after 
solving Eq. (U) for w and calculating 8w, which depends 
strongly on a competition between the different terms of 
Eq. (UJ at the z interval considered. For the sake of 



completeness, we also performed the same analysis for 
the other two cases (with 15% and 3% accuracy in H{z)) 
discussed earlier. In both, the scattering of the points is 
large enough to hamper any definitive conclusion on the 
expected evolution of the cosmic equation of state. 



V. CONCLUSIONS 

Cosmography explores the possibility of extracting the 
maximum amount of information from cosmological mea- 
surements, as well as the assumption that the universe 
can be modeled by the Friedmann- Robertson- Walker line 
element without assuming any dynamical theory to de- 
scribe it. In this paper, differently from many analy- 
ses that study the current phase of cosmic evolution in 
the context of a given cosmological model, we have in- 
vestigated the prospects for a cosmographic mapping of 
cosmic acceleration using two different projections of age 
determinations of passively evolving galaxies recently dis- 
cussed in the literature 0, [ljj. We have found that a 
model-independent check of the reality of cosmic accel- 
eration requires very accurated measurements of the ex- 
pasion rate (~ 1%), which seems to be far from the per- 
spectives of current planned surveys. As a step further 
in this analysis, we have also used the direct relation be- 
tween q(z) and w(z) to study possible constraints on the 
dark energy EoS from future observations (Fig. 4). 

Finally, it is also worth emphasizing that the method 
discussed here can be applied regardless of the technique 
used to obtain H(z) measurements. In this regard, we 
note that H(z) measurements from radial BAO methods 
(despite the systematic errors introduced mainly by dis- 
tortion effects) surpasses the age method in precision and 
can extend the H(z) measurements into deeper redshift 
ranges [f| (see also [l3| for expected measurements of ex- 
pansion rate at z < 2.5 by observing the Lyman- forest 
absorption spectra of high-z quasars). A detailed anal- 
ysis involving expected high-z estimates of H(z) from 
current planned BAO surveys is currently under investi- 
gation and will appear in a forthcoming communication. 
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